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ABSTRACT 

We investigate the screening effects on the nucleon-nucleon elastic cross section inside nuclear 
matter at zero and finite temperature. The N-N interaction is described phenomenologically 
via meson exchanges (a, it and uj) in the framework of a relativistic lagrangian model. The 
expressions for the in-medium meson propagators, which take into account the effects of matter 
polarization and renormalized vacuum polarization, are used in the calculation of the various 
meson contributions to the total cross section in the one-boson exchange approximation. The 
final expressions allow a satisfactory fit of the Coulomb substracted proton-proton cross section 
data in vacuum. At finite density or temperature the cross section is reduced as compared to the 
vacuum values and compatible with the existing experimental constraints. At high density the 
p-p differential cross section gets increasingly forward and backward peaked. The a-u mixing 
plays a key role in the understanding of the results. 
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(1) INTRODUCTION 

In plasma physics the screening effects of the medium modify the effective interparticle interac- 
tions which often become very different than in vacuum. The knowledge of the effective interac- 
tion is necessary for a proper analysis of the collective behavior of the plasma. In the particular 
example of nuclear matter, the screening of the medium on the N-N interaction (described 
in vacuum by a two-nucleon potential which has a Yukawa- like expression), leads to a new 
form of the effective potential which is long-ranged and oscillatory pi . Such " Friedel-like" 
oscillations can induce, under particular conditions, the appearance of some new observable 
properties in the nuclear plasma, which should be the analogues of the ones encountered in 
degenerate electromagnetic plasmas. 

Another important parameter related to the nuclear interaction, which can be modified by 
the screening, is the N-N cross section. It is a basic ingredient in the dynamical analysis of 
heavy ion collisions. Most of the early calculations in this field have been performed using a 
constant value of the cross section (a = 40 mb) or a fit of the experimental data in vacuum 
[||. However, there is an ever growing bulk of evidence that the modifications introduced 
by in-medium effects on the cross section have a non-negligible influence on the findings of 
these calculations. For example, a criterion which is often used to determine the stiffness of the 
equation of state is the amount of transverse flow. But the collective flow, and its disappearance 
at the balance energy, is also used to set constraints on the magnitude of the N-N cross section 
in the medium. It is therefore not possible to disentangle completely the effects coming from 
the equation of state or from a reduced cross section. It has lately been found by several 
authors, from BUU g §, |, 0, |] as well as QMD @ or AMD go| calculations, that a lower 
cross section in the medium is favoured by a detailed analysis of the transverse flow. Moreover, 
the cross section is an essential parameter in the calculation of transport coefficients and the 
knowledge of its screened behavior allows the determination of these coefficients as functions 
of the thermodynamical state. 

There exists several calculations [jllj - |19| | of the in-medium elastic N-N cross section. Only 
one of these calculations considers finite temperature [|l7]] as well as density. Most of them were 
done in the Brueckner formalism and are appropriate for the low energy range Ei^ ~ [50 — 300] 
MeV. Nevertheless, the relativistic effects should become important not only for higher energies, 
but (due to the medium effects) for high densities and temperatures. The purpose of this paper 
is to perform a fully relativistic calculation of the screening effects of the medium on the elastic 
N-N cross section in nuclear matter at finite density and temperature. 

An acceptable (phenomenological) description of the nuclear interaction in terms of meson 
exchanges requires the analysis of several (scalar, pseudoscalar, vector, pseudovector and tensor) 
meson couplings Efl|. This work is a step in a program of analysis of the screening effects on 
the complete nuclear interaction in nuclear matter, where all these couplings will be considered. 
Here we treat only the case of the exchange of scalar (a), pseudoscalar (tt) and vector (w) meson 
couplings. Although the present treatment of the interaction is not complete, it allows for a good 
fit of the experimental data of the elastic N-N cross section Q in vacuum and gives a satisfactory 
description of the screening effects on this parameter. Indeed, preliminary calculations including 
the Yukawa and derivative tensor couplings to the p meson, which will be published elsewhere, 

1 Here, we restrict to the Coulomb substracted proton-proton p-p or neutron-neutron n-n cross section 
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indicate that the corrections to the present calculations introduced by these suplementary 
couplings are small. Nevertheless, due to the isospin asymmetry associated to the p-N coupling, 
the proton-neutron (p-n) elastic cross section is more affected by the p-exchange contributions 
than the p-p and n-n ones. 

At highly relativistic energies, the A(1232) resonance and eventually more massive resonances 
such as A* (1440) are expected to play an important role. We calculated here the elastic channel 
of the N-N cross section and did not consider these resonances. The way in which the taking 
into consideration of the A resonance might modify our analysis will be discussed in the last 
section. Basically two types of corrections will have to be considered: 

First, the inelastic channels N + N -> N + A and N + N -> A + A open at T lab > 350 MeV. 
At energies higher than 800 MeV, the inelastic contribution represents 50% of the total N-N 
cross section. The A production and absorption via the inverse processes N + A —* N + N 
... etc mechanism is the input for the calculation of the subsequent pion production through 
A decay (A — > N + ir). The theoretical pion distribution and spectra obtained in this way 
serve as a basis for comparison with the values observed experimentally in relativistic heavy 
ion collisions. This important topic deserves a separate study and will be treated elsewhere. 

Secondly, the A's can also appear in internal loops. They will thus modify the screening of 
the interaction. There exists an abundant litterature on the modification of the ir dispersion 



relation by the A (see e.g. [21], [2^, gj| ). This topic is not yet completely settled, especially 
concerning the relativistic treatment. We argue that the screening of the pion will not play an 
important role for the proton-proton cross section calculated here, and will discuss the case of 
the modifications to the a and u dispersion relations in the last section (see section (6)). As 
far as we know, the modification of other mesons than the ir by delta-hole loops has not yet 
been considered in the litterature; this will be the subject of future work. 

In section (2) we give the lagrangian model, define the Hartree equilibrium and introduce the in- 
medium meson propagators at finite temperature, obtained in the one-loop approximation and 
including the vacuum polarization effects. Section (3) is devoted to the calculation of the elastic 
N-N cross sections from the screened propagators in the one-boson exchange approximation. In 
section (4) we determine the model parameters in order to fit the available experimental data 
on cross section in vacuum. In section (5) we present the results with the analysis of the effects 
of density and temperature on this variable, and conclude in section (6) with a discussion on 
the possible extensions of the model. 



(2) THE MODEL AND THE MESON DISPERSION RELATIONS 
2.1 Lagrangian density 

The cross section will be calculated by considering that the nucleons are interacting through 
the exchange of a, u and ir mesons. The dynamics of the model is given by the Lagrangian 
density 



L = 2 [ V* 7 " (^VO — (9^) • 7^ ] — m V'V' 
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+ -[{d»e){d u o)- ^ 2 a 2 \ 

- -jrF^ v - Itf&'&v 

+ ga^cf^ + gwipi' 1 tiiiip + ip j 5 TIT 1p (1) 

where 

F^ u = d»Qj v - d v u» (2) 
and the hat symbol over characters denotes quantum operators. 

A meson-meson interaction term of the form a.7? 2 is usually added to this lagrangian model in 
order to improve the fit of the pion-nucleon scattering amplitudes in vacuum. If the value of the 
associated coupling constant is properly constrained, cancellations of large contributions to the 
s-wave scattering lengths occur and the predictions of the model become reasonable in vacuum 
[p4|]. Nevertheless, no such cancellations occur at finite density and temperature, unless this 
coupling constant be differently constrained in every thermodynamical state. Moreover, for 
some values of the coupling strength this term can lead to difficulties related to the tachyonic 
poles introduced in the pion propagator by the imaginary character of the effective mass of the 
quasi-pion in the Hartree approximation. Finally, as we have tested, this term has no important 
quantitative incidence on the calculated cross sections in this approximation. Consequently, we 
shall not consider this interaction here. 

In the calculations performed in this paper, the numerical values of the constants appearing 
in the Lagrangian (|l|) are fixed as follows: The meson and fermion masses are fixed to their 
"physical" values p a = 550 MeV, = 783 MeV, fj^ = 138 MeV and m = 939 MeV. 
The solution of the model in the relativistic Hartree approximation (first order approximation 
in our scheme) leads to the equation of state of the plasma p5j| . Then, we chose for the u> 
and a coupling constants the values which give a satisfactory fit of the saturation properties in 
this approximation. These are g 2 = 183.3 (fia/m) 2 and gl = 114.7 (^/m) 2 . With these 
values, saturation is attained at Pjq = 1.40 /m _1 with a binding energy Eb = —15.85 MeV. 
Nevertheless, in the higher order approximations leading to the calculation of the in- medium 
elastic cross sections we shall often use different values of the coupling constants in order to 
fit the experimental values of this parameter in the vacuum (see section 4). This is justified 
because in our cluster expansion scheme [^] the values of the parameters must be redefined at 
every order. 



2.1 Dispersion relations 

Beyond the Hartree approximation, we treat the in-medium N-N interaction as mediated by 
dressed mesons. The one-loop to — a meson propagators at T = have been obtained for 
the first time in Ref. [27]. The dispersion relations and propagators of the mesons at finite 
temperature used here were obtained from a linear kinetic analysis of perturbations around the 
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Hartree ground state, starting from the relativistic quantum B.B.G.K.Y. hierarchy |p8|, |29||. This 
cluster expansion differs from the usual loop expansion, but one can show that the first order 
of the B.B.G.K.Y expansion coincides, at zero temperature, with the one-loop approximation 
[||. Nevertheless, the higher order corrections in the kinetic expansion differ from the results 
of many-loop calculations. 

Renormalized vacuum polarization contributions are crucial in obtaining a physically reason- 
able behavior of the propagation modes [29]. A first analysis of the structure of the one- loop 
a — io propagators at zero temperature has been performed by Lim and Horowitz |3(| in the 
semiclassical approximation (no vacuum polarization effects). It has been partially extended 
to the finite temperature case with vacuum effects in and (completed with a study of the 
the zero-sound modes and a improved choice of the renormalization point) in Ref. p2| . 

The dispersion relation matrix as obtained in Ref. [29| has the form 



D(k) 



-k»k v + (k 2 - nDg^ 

-g 2 X{k) 





-k 2 + ^1- g 2 a U a (k) 







k 2 



lil + glll^k) _ 



(3) 



where the general expressions for the polarizations H(k) of the meson fields as functions of the 
thermodynamical state (density and temperature) and 4-momentum k^ = (to,q) are summa- 
rized in appendix A. The ir propagation is decoupled from the dynamics of the other mesons 
because of parity conservation. There is a mixing between the a and longitudinal-a; propa- 
gation modes. The dynamics of the mesons in the medium is affected by the real as well as 
the imaginary parts of the polarizations. As mentioned above, the analysis of the propagation 
modes associated to the real components of the polarizations has been performed already. Let 
us make some comments about the behavior of the imaginary polarizations (see also appendix 
A.2): 

a) Except for the mixing terms, the baryon and antibaryon contributions to the imaginary parts 
of the polarizations have opposite signs. The antibaryon contributions vanish in the T = 
limit and the baryon contributions reduce to the usual expression for the matter contribution to 
the imaginary polarizations. The vacuum contributions to the imaginary parts have the same 
form as in the T = case and are finite. 



b) For time-like modes, there are two regions in the first sector of the q 
the parabola defined by 



to plane separated by 



ui = ^q 2 + AM 2 (A = 1 - AM 2 /(uj 2 - q 2 ) = 0). (4) 

Here, M = m — g a < a > is the effective baryon mass. As in the T = case the imaginary 
parts vanish in the region between this parabola and the diagonal u = q (A < 0), and the 
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modes are undamped there. Over this parabola (A > 0) the modes are damped by the decay 
into particle-antiparticle pairs. 

c) For spacelike modes in the T = limit, the imaginary parts vanish outside the region 
defined by 



\/W " 1? + M 2 -^p 2 F + M 2 < uj < \l (pp + q) 2 + M 2 — + M 2 (5) 

(pF is the Fermi momentum) where the modes are damped by the decay into particle-hole 
pairs. This region arises because, at T = 0, the baryonic Fermi factor in Eq.(f37|) becomes 
a step function. At finite temperature the imaginary parts are finite for any space-like mode. 
Nevertheless, for moderate temperatures, the behavior of the Fermi factors makes the imaginary 
parts small outside the particle-hole region. 



(3) THE O.B.E. CROSS SECTION 

The in-medium cross section is calculated by describing the nucleon interaction through the 
exchange of single dressed a, to and 7r mesons. The generic O.B.E. diagram is given in Fig.(l). 
The fermion lines correspond to the nucleons in the medium which acquire an effective mass 
and momentum from the ambient mean field in thermodynamic equilibrium. The properties of 
the effective nucleons are calculated in the Hartree approximation with renormalized vacuum 
fluctuations taken into account |33| . The dressed propagator matrix of the mesons is obtained 
by the inversion of the dispersion relation matrix (||). 



3.1 Dressed meson propagators 

The expressions of the polarizations given in appendix A were calculated, without loss of 
generality, in a frame (Rl) where the background fluid is at rest and the momentum tranfer 
is k^ = (u>,0,0,q). In the following we will be working in the center of mass system of the 
collision (R2) in which the matter comoving system travels with the quadrivelocity u^ 1 . It will 
therefore be more convenient to put the dispersion relations and propagators under a compact 
covariant form. 

We define the auxiliary quadrivector rf and the projectors onto the longitudinal {A^ v ) and 
transverse (T^ v ) modes as follows: 

i. „, 

(6) 
(7) 
(8) 
0) 



if 


= u^- 


k 2 




= <r- 


k»k v 
k 2 




= sT- 


rfrf 




rj 2 




rfrf 
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so that the polarizations now read 



n&, = n xf ,7h 2 ] 1/2 



(10) 



in terms of the longitudinal, transverse and mixing parts H^l, ^-uT and n x . These scalar 
quantities are related to the components IT^ , II* , 11°^ defined in appendix A by 



U. 



u 



-n 



00 . TT TT 

■u i lL u>T — 11 



11 . tt2 

ui i xx x 



(n; 



\2 
aui) 



The expression for meson propagator matrix becomes 







G» c 
V 



G a 
G n 



where the different elements (propagators) are given by 



(12) 



G a 



Pi + sl^L 



' (k 2 -fil+ g 2 U a )(k 2 - m£ + cfcSluL) + g^g^l 



(13) 



G^ 

w (TUI 



rj-ia> 



1 



k^k v 



k 2 - fJ% + g%Uu>T vlk 2 

k 2 -n 2 a + g 2 a U a 



(k 2 + g 2 U a )(k 2 - m£ + g 2 Ji^ L ) + gagjl 
k^k 



2 U2 



-IT}' 



g<jgu n> 



(k 2 - n% + g 2 a u a )(k 2 -1**, + g 2 Jiu L ) + g*gu n 2 x 

irfG* 



(14) 

(15) 
(16) 



G-jt 



1 



'k 2 -fil + 5 2 n ^ 



(17) 



In inverting the dispersion matrix Eq.(Q), the meson branches associated to the zeroes of the 
determinant become poles of the meson propagators. We are interested here in the spacelike 
region of the q — u plane and, when vacuum polarization is taken into account, there are only 
three poles of the propagator in this region, associated to the "tachyonic" branches of the 
mixed ^longitudinal, transverse and pion modes [29]. These branches come from the vacuum 
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polarization terms and arise at large k, where the one-loop calculations fails J34J. In fact, for 
these values of k (short distances) the point particle approach becomes inadequate and the 
structural properties of nucleons must be taken into account. In a more fundamental theory 
of pointlike particles, vertex corrections could be calculated self-consistently (see e.g. [j35[ ). In 
nuclear matter models, such corrections have been introduced in analogy with the QED case, 
with the suppression of high momentum transfer arising from the virtual Bremsstrahlung of 
neutral vector u> mesons |36|, 37]. In our phenomenological model, this can as well be done 
through the introduction of phenomenological monopolar form factors for the nucleon at every 
vertex in the boson exchange diagrams as well as at every vertex in the loop calculation of the 
boson propagators. This leads to new expressions which, in /c-space, are obtained from the 
original ones by multiplying every coupling constant by the corresponding form factor. These 
form factors have the generic form 



F(k) = (A 2 - ^ 2 )/(A 2 - k 2 ) (18) 

where fi is the scalar, vector or pion meson mass in every case and A is a cut-off parameter. 
After the introduction of the form factors, if the cut-offs are in appropriate ranges, the spurious 
"tachyonic branches" as well as the associated poles of the meson propagators on the uj < q 
region disappear. The renormalized vacuum contributions are crucial in obtaining physical 
(stable and causal) propagation modes p9| ]. If they are neglected, new unphysical spacelike 
modes are present at small q which can not be eliminated by the form factors and adulterate 
the pole structure of the propagators. 

3.2 Transition matrix and N-N cross section 

The O.B.E. amplitude A4 is obtained in the first Born approximation from the diagram of 
Fig.(l). In condensed form, we have 



M= (U 3 T m U 1 )G^ l (U A T n U 2 )- exchange (19) 

m,n=cr,ui,7r 

In this expression the U's are the spinors defining the incoming (U±, U2) and outgoing (U3, U4) 
nucleon states. The indices i = 1, 2, 3, 4 stand for the effective momenta (pf )* = pf — g^ < uj^ > 
and the spin-isospin quantum numbers S{, U. G™™ is the propagator matrix (|l^) which depends 
on the momentum transfer k^ = Pi* — P3* ■ T is the coupling matrix defined by 



(20) 



In the explicit calculation of the transition matrix from these formulae, all the coupling con- 
stants in G mn and T m are to be multiplied by the corresponding form factors. The exchange 
term is obtained by interverting the indices 3 and 4. 
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In the referential (Rl) where the fluid is at rest, u M = (1,0) and 

p» = P» + k»/2 ; p^ = K^-k^/2 ■ p£ = - 1^/2 ; P 1 = K^ + k^/2 



with 



^ = (^,0,0,9) 

In the center of mass of the collision (referential (R2)), the background fluid moves with a 
velocity = (7, yv) parametrized by 

' v sin a sin <p ' 

v= I v sin a cos 99 ] (21) 
v cos a 

The momenta of the incident and scattered particles are (energy conservation is already taken 
into account) 

p1 = (E,p) , p% = (E,-p) , = (E,f) , p» = {E-p') , 

the scattering angle is 9 defined by p.p 1 ' j \f(p 2 p 2 ) = cos#. The momentum transfers for the 
direct (£31) and exchange (k^x) diagrams are 



Ki = (o,P+f) 

One goes from (Rl) to (R2) by performing a boost of four-velocity (7,7*?) and a rotation. In 
the referential (Rl), 

k 3i( R1 ) = ( w e,fe) 
&41CRI) = (w©,<f©) 

a^/ e ^ = 7fp[±sinasin(^sin^ — cosa(l =F cos 0)1 

g f9) = [^ e ^ + 2p 2 (lTcos0)] 1 / 2 

The expression of the differential cross section for elastic N-N scattering in the center of mass 
of the two colliding particles is given by 



da/dn = 1/(64tt 2 s)\M\' 



(22) 



where s is the Mandelstam variable s = AE 2 . In the following we consider only the spin-averaged 
cross section obtained by replacing \M\ 2 by \M\ 2 = (1/4) S S1 S2 S3 S4 \M\ 2 in Eq. (|22|) . We 
obtain the p-p cross-section for the isospin quantum numbers t\ = t% = £3 = = 1 and the n-n 
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cross section for ^1 = ^2 = ^3 = ^4 = 0. a pp and a nn are equal in our approximation where the 
Coulomb interaction was substracted and the neutron and proton masses taken to be equal. 

The explicit formula for da/dQ is given in Appendix B in terms of the Mandelstam variables. 

In the next section the values of the coupling constants and the cut-offs of the form factors will 
be fixed in order to fit the cross section experimental data in vacuum. 



(4) CHOICE OF MODEL PARAMETERS 



4.1 Vacuum effects 

We recall that the vacuum effects were renormalized (see appendix A) in such a way that the 
contribution of the vacuum polarizations at vanishing temperature and density cancel on the 
mass shell of the mesons: 

U v i ac (k = hi) = at fe? = nf, i = {a, u, vr} (23) 

For other values of k however, the polarisation of the vacuum will bring a contribution to the 
free propagator 

G°dk) = - 2 1 , nUk) = - 9 mnr c (k) 

A 2 2 

with 9i (k) = ft-^ff (24) 

The taking into account of the vacuum polarization has the effect of reducing the effective 
k— dependent meson masses fj,i*(k) at high momentum transfer k. As mentioned above, the 
vacuum polarization contributions are essential in obtaining a reasonable pole structure of the 
propagator |2{|). 



4.2 Fitting procedure 

Our model includes only tt, a, u meson exchange. It is a minimal model which allows to 
reproduce the general features of the nucleon-nucleon cross section with only five parameters: 
the couplings of the scalar and vector mesons g a , and the cutoffs A CT , A w , A^. We fix the 
pion-nucleon coupling constant from experimental data on pion-nucleon scattering. Recent 
work (see e.g. |38|) quote a value slightly lower (g^/(47r) = 13.7) than the standard textbook 
value (<7^/(47r) = 14.4). We chose the latter one in order to be consistent with previous work. 
This will have no important incidence on our results anyway. The meson masses are kept 
fixed at their physical value [i^ = 138 MeV and [i^, = 783 MeV. The scalar meson is not 
physical but is a phenomenological way of describing correlated two pion exchange, so that the 
constraint on its mass is less important 0. We take the standard Bonn potential value ii a = 550 
MeV. The value of the iv-N-N coupling constant is not very well known experimentally. This 

2 The a might be identified with the fo (400-1200) broad resonance mentioned in the last issue of the Particle Data Book fji^ . 
The controversy about the identification of this scalar particle as a broad tt — n resonance was revived by recent observations 
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situation may be improved by future planned experiments near the u production threshold. 



Quark models |4l], 42|, [43|] predict g u ~ 8.5 — 8.8. From the decay width T(oj — > e + e~) and in 
the vector dominance approximation, a value ~ 8.6 can be estimated. On the other hand, 
the Bonn model value <?^,/(47r) ~ 20., obtained by fitting the nucleon scattering phase shifts 
is much larger, but it can be reduced by taking into account correlated tt — p exchange. We 
will let vary and will see in the next subsection that our fit yields values in the range 8-9. 
The cutoffs should be chosen larger than the mass of the meson exchanged. There is also an 
upper limit set by the constraint that the tachyonic modes of the dispersion relations should be 
removed (see the discussion of section (3.1) and Ref. (^9|). In practice, we restrict the search 
to the range 1.5 m < A» < 1700 MeV. 

Publication of experimental data is scattered throughout the litterature. Recent compilations 



and smooth fitting of the available data are to be found in Refs. [44, |4jj]. These compilations 
do not substract the Coulomb scattering in p-p data at low momentum exchange. Coulomb 
scattering is most important at low energy Ti a b < 40 MeV and/or scattering angles < 10° — 
20°. We consider here a model of the nuclear interaction only, since the cross section we are 
calculating is intended for use in the current models of nuclear matter (BUU, QMD . . . ), where 
the Coulomb interaction can be neglected. Simple fits of Coulomb-substracted experimental 



data are given by Cugnon et al. 0, 46] for use in their cascade model. Some Coulomb-subtracted 



data points are also given in the older paper of Chen et al. [47]. 

We performed a simultaneous fit of the total elastic cross-section and the differential cross 
section at 9 distinct values of the cm. energy. The fit was obtained using a weighed least- 
squares method where the data points are chosen as follows: For the total cross section a pp , we 
used 12 points distributed over the energy range £ [40—1000] MeV. For the differential cross 
section, da pp /dQ, we used 8 points distributed in the range of scattering angles 9 G [10° — 90°] 
for each value of the energy T lab =60, 100, 250, 350, 460, 560, 660, 800, 1000 MeV. Fitting 
only dapp/dQ would in principle be enough since a pp is obtained by integrating over 9, but we 
obtain better results by taking a pp -data with a weight of 10 % - 20 % into account. Very-low 
or very-high energy data for da pp /dQ was also attributed a lower weight. We found in this way 
X 2 values in the range 10 - 12. 

Due to the reduced number of free parameters, our model is not be expected to do as accurate 



a job as the last Bonn |20|, |48|, 149(1 , Nijmegen |50|] or Argonne |)1|] potentials for example. 
Nevertheless it will be seen in the remaining that a reasonably good fit can be obtained for the 
p-p (and, which is equivalent here, n-n) total and differential elastic cross section. It is more 
difficult to reproduce the p-n cross sections. As a matter of fact, the T = 1 isospin channel 
which is necessary for the calculation of p-n scattering requires the presence of the p-meson 
exchange term. For this reason, in this paper we show only results for the p-p cross-section. 
Inclusion of more mesons (p, n, 5) and calculation of the p-n cross section is currently under 
way. 

4.3 Fitting parameters 

Depending on the details of the fitting procedure, we found several equivalent parameter sets. 

The total elastic p—p cross section a pp is well reproduced in the energy range [40 — 1000] MeV by 
the parameter sets (A,B,C). (see Fig. 2). We also quote the somewhat extreme choices (D,E,F). 
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9a 


9w 


A CT 


A w 


A^ 


set A 


3.80 


9.31 


1298.8 


1240.5 


362.1 


set B 


3.27 


8.38 


1622.0 


1324.7 


355.2 


set C 


4.87 


8.99 


951.3 


1255.3 


397.4 


set D 


7.93 


8.93 


703.3 


1237.3 


376.1 


set E 


8.9 


37.7 


1643.2 


903.2 


587.2 




m a 


9uj 


A CT 


A w 


A^ 


set F 


763.3 


8.06 


912.2 


1348.5 


390.2 



By comparing our results at finite density and temperature, this will allow to investigate the 
sensitivity of our conclusions to the choice of model parameters. (D) is obtained by fixing the 
coupling constants to their saturation values (see section (2)), so that it contains only the three 
cutoffs as free parameters. (F) is obtained by fixing the coupling constant of the a meson 
equal to that of the pion and instead let vary the a mass. Sets (A,B,C,D,F) all correspond 
to the same minimum in parameter space. Sets (D) and (F) have to pay for a too high g a 
by a correspondingly low A CT . Set (E) corresponds to another (not very stable) minimum in 
parameter space. It is somewhat better at low energy but grows rapidly bad at energies higher 
than 500 MeV. There a very strong vector coupling is castigated by a small A w . The best 
choice is set (A), then we have (B) > (C) > (D) > (F) > (E). In all cases we notice that the 
cut-off parameter of the pion is smaller than usually admitted in Bonn parametrizations. This 
is due to the absence of compensation by the p meson. Preliminary results indeed show that 
the inclusion of the p meson allows higher values of A,,-. Still the A^'s of all sets shown here are 
always above twice and a half the pion mass. 

The differential cross section is somewhat more difficult to fit accurately (see Fig. 3). It is 
possible to perform a good quantitative fit at high energies. At lower energies, the fit becomes 
qualitative. The same remarks about the goodness (badness) of the parameter sets as for the 
total cross section apply here. The low energy behavior would be improved by increasing g a . 

If we look at the contributions of each meson separately, we see that the value of the cross 
section is essentially driven by the a and uj mesons. In a way similar to what happens for the 
nuclear potential p3], the value of a pp results from a delicate balance of a and u contributions. 
The pion describes the low momentum-transfer, long-range part and the detail of the differential 
cross section at forward angles. 



(5) RESULTS AND DISCUSSION 

As a first important statement, it is found that the various choices of fitting parameters essen- 
tially give the same behavior at finite density and temperature. We display here results for the 
set of parameters (A), and comment about the other sets when appropriate. 

5.1 In- medium cross section from screening 

The total elastic p-p cross section at zero temperature is displayed in Fig. 4 for the parameter 
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set (A) and for various values of the density. It is seen that it decreases with increasing density. 
The reduction is larger at smaller momentum transfer and energy, so that the total cross section 
becomes approximately constant throughout the whole energy interval for n ~ n sat . Further 
increasing the density does not much modify a pp . 

We do not observe significant differences for the total p-p cross section calculated with param- 
eters sets (A,B,C,D,F). For sets (A,B,C), o pp (n > 2 n sat ) settles around ~ 17 mb for all T^. 
For sets (D,F), we have a somewhat lower value a pp (n > 2 n sa t) — 15.5 mb, but we saw that 
both models had a too high scalar coupling g a and a too low scalar A CT cutoff. The general 
features (low-density reduction, flattening to a T^-independent value, constant value at high 
density) are also obtained with set (E), but it settles to a higher value a pp (n > 2 n sat ) ~ 71 
mb. This high value is to be traced to the unreasonably strong g w and low cutoff K u . 

We also calculated the total p-p cross section at finite temperature. The ratio of the value of 
the cross section in the medium to its value in the vacuum is shown in Fig. 5 as a function of 
density n/n sat for seven values of the temperature T=0, 40, 80, 120, 160, 200 and 240 MeV and 
for two values of the energy of the beam energy T} a f, = 100 MeV and l} a b = 300 MeV. This same 
ratio is displayed as a function of temperature on Fig. 6a for five values of the density n/n sa t 
= 0, 0.1, 1., 2. and 4. for the same values of the beam energy. These results are summarized 
in Fig. 7 as a contour plot in the density-temperature plane for incident particle energies 
= 100, 300, 500 MeV. 

At n = 0, a similar behavior is obtained with increasing the temperature as at T = with 
increasing density (the a pp cross section decreases at high T) but the temperature dependence 
is not as strong as the density dependence: There is almost no difference between the T = 
and T = 120 curves. The same quantitative effect is obtained by heating the vacuum at T = 
(150, 170, 200) MeV, or by compressing cold matter at n = (0.1, 0.5, 1) times the saturation 
density. The cross section at n = becomes flat over the energy range [40 — 1000] MeV at 
T=250 MeV. 

At fixed non vanishing temperature, the cross section is a nonmonotonous function of the 
density as can be seen on Fig. 5. This figure also reveals that the cross section is more sensitive 
to temperature in the density range of interest for nuclear collisions n/n sa t= 0.2 - 2. 

At fixed non-vanishing density, the cross section is non monotonous as a function of the tem- 
perature. It first slightly increases, reaches a maximum and decreases rather rapidly at high T. 
The turnpoint temperature increases with increasing density. In brief, the behavior of the cross 
section a pp with varying T qualitatively follows that of the nucleon effective mass (compare 
Figs. 6a and 6b). 

The remarks done on the example of parameter set (A) as to the behavior of the total cross 
section with temperature hold for (B,C,D,F) parameter sets as well, save some unsignificant 
quantitative differences. 

Let us finally discuss the density and temperature dependence of the differential p-p cross section 
(see Fig. 8 and Fig 3c). Since da pp /dfl is quantitatively not so well reproduced as the total one, 
we have to be more cautious in our conclusions. We nevertheless notice some general features: 
at high density or temperature, the differential p-p cross section becomes very forward peaked 
(or equivalently backward peaked since the p-p cross section is symmetric around 9 = 90°). 
This trend is not affected by the same "saturation" effect as for the integrated cross section. 
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At low density on the other hand, forward scattering is slightly suppressed with respect to the 
free value. 

The behavior of the cross section as a function of density and temperature can directly be 
related to that of the polarizations. As a matter of fact, the polarizations appreciably differ 
from zero when the momentum transferred by the meson is lower than the threshold value 
\/—k 2 = 2 pp for particle- hole production/recombination ^. When this occurs, the mixing 
between the sigma meson and the longitudinal mode of the omega meson through particle-hole 
loops is important, and the balance between the attractive and repulsive parts of the interaction 
is disturbed. 

Let us recast the propagators of the mesons inside the medium in a form as close as possible 
to the free one (e.g. G° a = l/{k 2 - [/x™ c (A:)] 2 )): 

K ■ 

k?-\ M ,f(k)]2 + iiJ,f(k)T a (k) { 1 

K ■ 1 
Q = mix \v j T^ v 

k 2 ~ + irf L (k)T uL (k) + k 2 - [nf T {k)f + inf T (k)T„ T {k) K > 

\pf{k)f = £-9lKeIl a (n,T,k) (27) 
1AV = tif/gllmllvin^k) (28) 

(and similarly for ir, loL and uT), and 

K (k 2 + 9l^ L ){k 2 -/£ + g 2 a U a ) 

(k 2 - l , 2 +g 2 U ujL )(k 2 -^+g 2 U a )+g a g uj n 2 x [ 1 



with 



This defines momentum dependent "effective masses" and widths, and a mixing parameter. 
The ratio of the "effective masses" to their vacuum values and the mixing parameter are shown 
on Fig. 9 for T = and n/n sa t = 0.1. Also shown on this figure are the dynamical masses, 
defined as the solution of the dispersion relation at zero 3- momentum k = (ujf l ,0): fif yn = u>f , 
D(k) = 0. It can be seen that the mixing parameter deviates strongly from its vacuum value 
( = 1) as long as the momentum transfer is lower than 2 pp. The transverse part of the u 
propagator remains almost unchanged. On the other hand, the "effective masses" of the a and 
longitudinal part of the u follow the behavior of -Kmix> the effective a mass being reduced and 
the effective longitudinal lu mass enhanced in the same proportion. 

The mixing parameter comes as a multiplicative factor in G a and G^l , so that the mixing could 
be absorbed in the redefinition of effective couplings for the a and the longitudinal u, and be 
viewed as a weakening of the strength of these effective interactions. 

Finally, the non-vanishing imaginary part of the polarizations which describes damping of the 
modes through particle-hole decay, contributes a positive term in the denominator of the square 
of the propagators and further reduces the value of the cross section in the medium. 

3 The condition \/—k 2 = 2 pf is the same as that leading to the Kohn singularity at the threshold for particle-hole production 
which was at the origin of the Friedel oscillations in the N-N potential ll| . 
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It is now easy to understand the behavior of the differential cross-section at finite density. Since 
the momentum transfer is given in terms of the incident energy E = Ti a b + M and scattering 
angle 9 by k 2 = —2 (E 2 — M 2 )(l ± cos 9), the condition \/—k 2 < 2 pp yields at constant energy 
the 9 range affected by medium effects: 



cos^|>l-^L ( , )()) 



This corresponds to the range where da pp /d£l becomes very forward peaked (and also backwards- 
peaked, since da pp /d£l is symmetric by 9 — > ir — 9). 

5.2 Comparison with Brueckner calculations 

The calculation we performed here amounts to consider loop contributions, while the Brueckner 



calculations of Li and Machleidt [16], Faessler et al. [11], Aim et al. [17] consider ladder 



contributions. In fact, both types of contributions should be included. This is however still 
somewhat beyond the present state of the art, as a naive ladder-resummation of our screened 
propagators would lead to inconsistencies (see e.g. a discussion in |]52| ). 

Still it is interesting to study separately the effect of loops or ladders. The Brueckner resum- 
mation is expected to be more appropriate for low densities, since it is based on an expansion 
in na 3 (a stands for the range of the interaction). On the other hand, our loop resummation 
is expected to be better at higher densities []53| and energies. In the case of the Brueckner 
resummation, the corrections due to finite density and temperature come in as Fermi factors 
on the intermediate (real) nucleon states on the uprights of the ladders, while in the case of 
loop resummations, the Fermi factors are applied on the intermediate (virtual) nucleon states 
forming the particle-hole loop. Both models, although they consider very different physical 
effects, lead to very similar results: decreasing a pp with increasing density and saturation for 
n > 2 n sa t, flattening over the energy range. 

There are also some differences: we observe a modification of the in-medium cross section at 
high energies too, although the effect is smaller than at low energy. We also have a more 
rapid decrease with density for low values of n/n sa t- This result should probably not be taken 
too seriously: First, because, as already pointed out above, our model is not optimized for 
this range, and second, because these densities at T = fall within the domain of spinodal 
dP/dp < instability. Indeed, such a rapid variation of the cross section with the density is 
not observed when the background matter is heated moderately (T ~ 20 MeV) so as to come 
out of the unphysical spinodal region. 



5.3 Comparison with experimental constraints 

During the last five years, it has become a general trend to investigate the influence of a modified 
N-N cross section in the medium in numerical simulations in heavy ion collisions. We mention 
here, as a non-exhaustive sample, the work of Klakow et al. [I| or Aim et al. || on BUU codes, 
of Khoa et al. ||] on a QMD code, or of Tanaka et al. [54] with the antisymmetrized AMD 
code. 
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All these determinations are strongly model-dependent, but the fact that they all predict 
roughly the same effect, i.e. a 20 % reduction of the in-medium cross section at n = n sa t, 
although they rely on a very different physical analysis of the processes involved, sounds very 
encouraging. 

We compare in Fig. 10 our results as well as the Brueckner results of Li and Machleidt pi] to 



the BUU prediction of Klakow et al. Q as well as the AMD prediction of Tanaka et al. [54]. 

Tanaka et al. effected simulations of p+nucleus collisions using their antisymmetrized version 
of quantum molecular dynamics (AMD), and extracted a reduction factor for the in- medium 
N—N cross section with respect to its free value by comparing the results of the simulation to 
available experimental data in the energy range [0 - 200] MeV. The p+ 9 Be collision probes low 
densities in the surface, while the p+ 12 C, p+ 27 Al, p+ 40 Ca, which all give similar results, are 
thought to probe the nuclear density n = n sa t. The result of Tanaka et al. is represented on 
Fig. 10 by triangles for p + Ca and diamonds for p + 9 Be. 

Klakow et al. estimated the medium-modified N — N cross section by performing simulations of 
heavy-ion collisions A + A with a Boltzmann-Uehling-Uhlenbeck transport model. The heavy 
ion collision takes place with an initial energy per nucleon 200 MeV and the system is evolved 
until freeze-out. The information on the cross section is obtained from an analysis of the 
dependence of nuclear transverse collective flow on beam energy, and the measurement of the 
so-called balance energy Ef, a i at which the flow disappears. This energy indicates a transition 
in which one goes from attractive to repulsive scattering. The authors assumed a density 
dependence o m edl & free = 1 — Oi(n/n sat ) i.e. by taking the first order of a Taylor expansion at 
constant beam energy and zero temperature. The best agreement with experimental data was 
found to correspond to a = 0.2. The result of Klakow et al. is shown on Fig. 10 by two arrows 
for densities n = 0.5 n sat and n = n sat . 



The Brueckner calculation of Li and Machleidt [16] was described in the preceding paragraph. 
We display in Fig. 10 their result for the p-p cross section at n = 0.5 n sa t (dotted line) and 
n = n sat (continuous line). 

Our results for the in-medium p-p cross section are represented on Fig. 10 by four dashed 
and dot-dashed lines at densities n = 0.1n sat and n = n sat and for temperatures T = and 
T = 50 MeV, as limiting values of what density and temperature conditions can be reached in 
the collisions considered by Klakow et al. or Tanaka et al.. 

We can see that there is a fair overall agreement between these four approaches. Experimentally 
as well as theoretically, the cross section is found to be reduced with respect to its free value. 
There is also a qualitative agreement in the behavior with the beam energy: as the energy 
increases, the reduction factor decreases and saturates to a constant value. Both theoretical 
approaches (Brueckner, ours) overpredict the reduction of the N-N cross section with respect 
to the experimental estimates. Two reasons may explain the observed discrepancies: 

(i) The theoretical approaches consider that the collision is happening inside infinite nuclear 
matter, while in the experimental setup and the numerical simulations finite size effects play 
an important role, so that in a dynamical situation the dressing of the mesons may not be fully 
effective. 

(ii) The calculation of Li and Machleidt as well as ours assume that the background matter is not 
far from local equilibrium, while the initial state of the collision is strongly out of equilibrium, 
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which also modifies the dressing of the mesons. 

It is interesting to note that the discrepancies observed with the T=0 calculations (Brueckner 
and ours) are reduced when finite temperatures are considered. This can be an indication that 
such finite temperature effects play a role in moderating the reduction of the cross section. 



(6) CONCLUSION AND OUTLOOK 



6.1 Concluding summary 

We have analyzed in this work some aspects of the screening effects on the nucleon-nucleon 
interaction, as due to the modifications introduced by the nuclear medium (at zero and finite 
temperature) on the propagation of the mesons which mediate this interaction. We applied here 
our results on dressed meson propagators to the calculation of the modification of proton-proton 
(or neutron-neutron) cross section for collisions occuring inside dense and hot nuclear matter. 
We found that the p-p or n-n total elastic cross sections were reduced with increasing density. 
The effect is more important at low energy, leading to a nearly constant cross section a pp ~ 17 
mb at all energies for densities above the saturation value. The temperature dependence was 
found to be less pronounced, with first a slight increase and then a reduction with respect to the 
free value. The differential cross section becomes increasingly forward (and backward) peaked 
as the density or temperature grows. 

The general features of the behavior with density are found to be the same as from a calculation 
with the Brueckner formalism effected by Li and Machleidt 16], although both approaches 
are very different. The reduction with density and its behavior with energy are compatible, 
although somewhat overestimated, with respect to the information extracted from BUU and 
AMD analysis of experimental data. Taking into account the effects of finite temperature in 
our approach seems to improve the agreement. 



6.2 Improved meson exchange picture 

While the calculation presented in this paper was restricted to the nucleon-nucleon scattering 
between identical particles through the exchange of tt, a and uj mesons, it is also possible to do 
the same for the neutron-proton cross section. In order to describe the T = 1 isospin channel, 
it is necessary to consider also the exchange of the p meson. This calculation is presently under 
way, as well as the inclusion of the r/ and 5 mesons of the Bonn OBE potential. 

Other interaction Lagrangians can be implemented using the same formalism, in particular with 
(non linear) meson interactions or chiral symmetry. A chiral symmetric interaction, although 
it is desirable from general considerations, has not yet reached the level of accuracy of the 
popular one-boson exchange picture in the description of A-A scattering, so that we did not 



consider it in this work. This field however is in rapid evolution |56| , 57, [5^]. For example, the 
Nijmegen group [58] quote a value of x 2 of 1.75 for a fit of the N-N scattering data with 12 
free parameters. It will bring a valuable piece of information to see how chirality modifies the 
behavior of the cross section, in particular in the high density or high temperature range where 
the chiral restoration takes place. 
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6.3 Background velocity 

Also under consideration is the case when the background nuclear matter is moving with a 
non-zero velocity with respect to the center of mass of the collision. In this case, the behavior 
of the propagator in other regions of the tu-q plane will be explored. For a given value of the 
energy of the colliding particle and of the scattering angle, the values of the energy to and 
impulsion q entering the expression of the polarizations is given in section (3.2) as a function of 
the relative velocity. In particular (see Fig. 11), the region of the zero sound branches can be 
crossed. Such branches arise from the mixed a - longitudinal u part of the dispersion relation 



, 32|. If (u),q) coincides with a zero-sound solution k£ s = (u) zs ,q zs ) of the dispersion relation 
k 2 -fi 2 a +glU a (k)) (k 2 - fJ 2 UJ + g 2 n ujL (k))+g a g uj n 2 < (k) = 0, (31) 



there will be a critical value of the relative velocity v cr u for which the real part of the dispersion 
relation vanishes. For relative velocities near v cr u the real part of the dispersion, which appears 
in the denominator of the propagators, is reduced. As a consequence, the cross section will be 
enhanced. 

A more detailed analysis of this and other effects of the relative velocity of the center of mass 
of the collision with respect to the background will be published separately. 



6.4 A resonance 

An important extension of the model should be the inclusion of the A resonance. The A 
and higher mass resonances will appear at several levels of the calculation: not only as the 
product of inelastic scattering processes, but also in a virtual A-hole pair which will modify 
the propagation properties of the mesons and in particular of the pion. At high densities and 
temperature, the nuclear matter turns into the so-called resonance matter with a softer equation 
of state. 

At the relativistic energies we are considering, a large number of A(1232) are expected to be 
produced through the inelastic processes N + N — > A + N and N + N — > A + A, and reabsorbed 
by the inverse reactions. In the vacuum, the inelastic channel begins to open for laboratory 
energies > 350 MeV; at > 800MeV, the elastic and inelastic contributions to the total 
nucleon-nucleon cross section are roughly of the same size. The production vs. absorption 
rates of A resonances is an important input in the calculation of the rapidity distributions 
dN^/dy and dNx/dy of the pions and kaons since these particles are produced through A 
decays A + A -> A + A -> iV + iV + 7r and N + N -> A + A -> N + A + K. These 
distributions are the quantities actually accessible by experiment and are used to characterize 
the thermodynamical state of the piece of excited nuclear matter formed in a relativistic heavy 
ion collision. Scatterings N + A —* A + A will also contribute to the transport properties. 

It is in principle straightforward to extend our formalism to calculate the N + N — ► A + A 
cross section once the dressed meson propagators are known. The A is coupled to the nucleon 
through A-N-7T and A-N-p vertices. The scattering matrix is obtained from terms such as 

Tr[G N (l)T vNN G N (3)r vNN ] Tr[G N (2)T nNA G A (4)r nNA ] ^(31)^(31) (32) 

... etc and the A decay width may in a first approximation be taken into account by integrating 
over a mass distribution function of the Breit-Wigner form. In this expression, Gn and G A 
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are the dressed, on-shell propagators of the nucleon and A respectively, and D n is the screened 



pion propagator. Similar calculations were performed already (see e.g. [58] in the vacuum 
and p0| , |6lU with dressed pion propagators). Nevertheless, before doing this, we think that it is 
profitable to perform a more thorough study of the consequences of screened meson propagation 
on the nucleon-nucleon interaction on the example of elastic scattering. 

An other way the presence of the A can influence our calculations is in the description of the 
background nuclear matter. The A can also appear in an internal line and contribute to the 
dressing of the mesons which mediate the interaction. It cannot be excluded that the corrections 
brought about by the creation of delta-hole loops in the propagators of the mesons might be 
sizable. But due to uncertainties and ambiguities in the description of the A, a prediction at 
the quantitative level seems premature. 

The quantum hadrodynamics Lagrangian of Serot and Walecka can be extended in order to 
include the A resonance. The resulting (nonrenormalizable) Lagrangian (see e.g. |6^]) intro- 
duces the A with the Rarita-Schwinger spinor The A interacts through the exchange of 
a and u like the nucleons, moreover the A couples to the nucleons through A-N-7T and A-N-/3 
vertices. A first uncertainty comes from the fact that the coupling constants of the A to the a, 
u and p are not very well known, and there is not enough data to constrain these parameters. 
A more serious difficulty arises from the Rarita-Schwinger formalism itself. As a matter of 
fact, the propagator ^ A ^ Af j, contains not only spin 3/2 but also spin 1/2 parts. There has 
been a controversy whether or not one should project out the 1/2 part and how. When the 
A is interacting, the theory becomes plagued by so-called off-shell ambiguities, i.e. new free 
parameters may be introduced at the vertices. Studies on the influence on the A propagation 



noticed a strong dependence in these off-shell parameters [22 



There exists an abundant litterature on the modification of the tt dispersion relation by the 
A (see e.g. [21, 22, |23| ). This topic is not yet completely settled, especially concerning 
the relativistic treatment as discussed above. We have seen that the pion did not play a 
very important role in the proton-proton cross section. As a matter of fact, we checked this 
statement by simply replacing the dressed pion propagator by the free one, and observed a 
modification amounting to 5% at most. For the neutron-proton cross section nevertheless, the 
pion propagator will contribute a sizable part, as well as the p. 

On the other hand, the modifications to the a and uj dispersion relations should be investigated 
in the case of pp scattering as well as for np scattering. Additional polarizations are given by 
terms such as: 

J d 4 p Tr[G A (p + k/2)T AAa G A (p - k/2)T AAa ] (33) 

As far as we know, the modification of other mesons than the pion by delta-hole loops has not 
yet been considered in the litterature; this will need to be investigated in future work. 



6.5 Form factors 

The phenomenological form factors account for the fact that the nucleons are not pointlike, 
but have a quark gluon substructure. As a consequence, very high momentum transfer q would 
probe very short ranges at which our meson exchange picture fails. The cutoff suppresses these 
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unphysical contributions. The standard monopole form chosen in this work (see Eq. flOf)) fits 
well the known data on the form factor in the vacuum. 

At finite density and temperature, we kept this form factor unchanged according to our level of 
approximation. In fact, the form factor is expected to be modified in the medium: the nucleons 
"swell" (EMC effect) with increasing density. In a simple bag model picture, this is due to a 
reduced bag pressure. Higher order calculations of the vertex function |57|, where the nucleon 
size is generated by a surrounding meson cloud, yield the same behavior. Higher momentum 
exchanges will thus be more strongly suppressed at higher density. 



APPENDIX A: Analytical formulae for the polarisations 

Without loss of generality, the explicit form of the polarization terms can be given in the 
particular frame where k = (o>, 0, 0, q). In this frame the only non- vanishing components of the 
polarizations are [29]: 



°(k) ; Hk) ■ n w 00 (k) ; u u 03 (k) ■ n w 33 (fc) ; 

IE, n (k) = n w 22 (k) ; IL a (k) ; U^k) (34) 
The expressions of these components are directly obtained from 



U a (k) = 8vr J d A p [(p-fc/2)(p + jfe/2) + M 2 ] .E(k,p) 
U au> ^(k) = WttM J d 4 pp^.E(k,p) 



(35) 
(36) 



U u ^(k) = 16vr/d 4 p [p^p v - k^k v /4 - (1/2) (p 2 - k 2 /4 - M 2 )g^ v ] .S(fe,p) (37) 



U w (k) = 8tt J d A p [(p-k/2)(p + k/2) - M 2 J .T,(k,p) 



(38) 



where 



E(k,p) ={f(p + k/2) - f(p-k/2)]/(k.p - ie) (39) 

and the ie term comes from the boundary prescription for the fermionic dressed propagator 
which enters the loop expressions of the polarizations. It gives the imaginary parts of these 
polarizations through the well know formula \/{x — ie) = V /x + iir5(x), where V is the principal 
part. The fermion distribution function f(p) in ([39]) is given by 
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f(p) = (2vr)- 3 5{p 2 - M 2 ). [fl + + Q 



H(-po)} 



(40) 



the f2's are the Fermi factors defined as : 



H(± P0 ).{1 + exp(/?[E(p) =F fi])} 



i 



(41) 



and H(x) is the Heaviside step function. [5 is the inverse temperature and /x is the auxiliar 
chemical potential, related to the true chemical potential by [itrue = H — 9w^o- The 
quasi-particles are on the mass-shell: 



A.l Real part of polarization tensors 

The matter contributions to the real part of the polarization tensors is inmediatly obtained 
from Eqs. (|35| ) to ( |38|) by excluding the last term in Eqs. (f40|). This last term corresponds to 
the quasi-nucleon distribution in the Dirac sea (vacuum distribution) and leads to divergences 
in the calculation of the real parts of the polarization tensors which must be renormalized. For 
the real part of the scalar polarization U a (k), the renormalized vacuum contribution has been 
obtained as usually, by regularizing the infinite vacuum polarization (which contains divergences 
proportionnal to k 2 , a and a 2 ) and adding to the Lagrangian the counterterms which allow 
to cancel such divergences. The finite undefined constants in the renormalized expression are 
determined by imposing a on-mass-shell criterion to the real renormalized dispersion relation 
for the scalar mesons at vanishing density |^9fl . The final form of the vacuum finite (v.f.) part 
of the scalar polarization is 



M = m — g a a 



(42) 



= (2k 2 - 8m 2 )6 a + 8(M 2 - k 2 /A)6(k 2 , M 2 ) 
+ 8(m 2 - jj/fj 2 /4)(^ CT 2 - k 2 )6 ak + (6M 2 - k 2 ) ln{M/m) 2 
2g 2 a 2 [6 + p a 2 /m 2 + 16m 2 ^ m + 8(m 2 - p a 2 /4)m 2 8 amm ] 
- (M 2 - m 2 ) [6 - /v 2 /m 2 + 89 a + 8(m 2 - ^ 2 /4)9 am ] 



(43) 



and the function 6(k 2 ,M 2 ) is given by the integral: 




(44) 



with : 



y = 1 - k 2 /(4M 2 ). 
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In writing down Eq. ([43|), we have introduced the notations: 



= e(fi a 2 ,m 2 ), 
e am = de(k 2 ,M 2 )/dM 2 , 
e ak = de(k 2 ,M 2 )/dk 2 

0„mm = d 2 9(k 2 ,M 2 )/d(M 2 ) 2 . (45) 

All these derivatives are to be calculated at the point k 2 = /U 2 , M = m. 

As it can easily be seen, there is no vacuum contribution to the mixing polarization part 

The finite vacuum contributions to the vector polarization tensor 11^ ^ u (k), is obtained from 



Eq.(37) using dimensional regularization and compensating the divergent contributions through 



the appropriate counterterms [29]. The final expression is 



U^f ^(k) = (4/3?r) [k 2 \n(M/m) - 2M 2 + (2M 2 + k 2 )6(k 2 , M 2 ) + Y] A^(fc) (46) 
where the projector A fMU (k) is defined by: 

A^(k) = - k»k u /k 2 (47) 

Here 9(k 2 ,M 2 ) is the same function defined in Eq. ( f44| ) and T is a finite constant which, 
determined through the on-shell renormalization conditions, becomes: 

V = 2m 2 - (2m 2 + ^ 2 )6{^ 2 ,m 2 ) (48) 

The renormalized vacuum contributions to the real part of the pion polarization obtained from 
the same on-shell prescription is [29 

IW(A;) = 2A; 2 [^ - 6(k 2 ,M 2 ) + ^ 2 9 wk ] - 2(i„ 4 0„ k 

+ (2M 2 - k 2 ) In (M/m) 2 - (2 - 2 /m 2 - 2/^ 2 7rm ) (M 2 - m 2 ) 

- 2g 2 a 2 (2 + ^ 2 /m 2 - 2^ 2 m 2 Q m ) (49) 

where the definition of the derivatives of the functions 9 is the same as in Eqs. ([W]) with the 
replacement of \x a by 

A. 2 Imaginary part of polarization tensors. 

The only non-vanishing components of the imaginary parts of the polarization tensors are: 
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n*,i(*0 ; (50) 

n .,/(fc) ; (5i) 

nL,/(*0 = (^)nL,/W; (52) 

DgKfc) = (w/^DgrCfc) = W^fc) ; (53) 

n&W = n£r(*) (54) 



the relations between the vector-mixing and tensor components being a consequence of the 
baryonic charge conservation. The expression for the finite temperature imaginary part of the 
scalar polarization is: 

roc 

Tl aI (u,q) = A(lo 2 - q 2 )/q / dE[H{A + ) - H(A~~)] [Q + - DT - 1] (55) 

JM 

where 



A ± = qVE 2 — M 2 — \loE±(uj 2 - q 2 )/2\ (56) 
0± = {1 + exp((3[ET fA)}' 1 (57) 

A = 1 -4M 2 /(w 2 - q 2 ) (58) 

The three terms inside brackets in Eq. ( pq) correspond respectively to the contributions of 
baryons, antibaryons and vacuum. Here and in the next formulae the vacuum contributions 
are finite. 

The imaginary part of the pion polarization is obtained from the imaginary part of the a 
polarization through 



Hn.i(u),q) = U a i(u,q)/A 



(59) 



The independent imaginary component of the mixing polarization vector is 



A ]\/f roo 

n2w>,?) = — dE 

q Jm 



E--)H(A- 



LO 



E+^)H(A^ 



(n + + or) (60) 



where the vacuum contributions vanish. 

The independent imaginary components of the to polarization tensor are 
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ttOO 

LL LU,I 




H{A~) 



E + - 
2 



H(A+) 



+ 1 



(61) 



and 



dE 



M 



M 



H(A-) 

h(a+)\ (n- 1 



0" + 1) (62) 



APPENDIX B: Analytical expression of the cross section 

The expression of the spin-averaged nucleon-nucleon cross section results from taking the square 
of the transition matrix Eq. fll9| ) and using the relation J2 S u (P-> s )^(p> s ) = 7-P + m - We use 
the short notation Si = 7-Pi* + m*. 

The calculation of the transition matrix is straightforward. We give here the result as a function 
of the Mandelstam variables and of the relative background velocity. Additional terms appear 
due to a-uj mixing, for example for the direct term: 



|M| 2 {^31*) -^31*) (Tr[5 37 ^ l7Q ] Tr[S 4 S 2 ] 

+ Tr[5 3 7 M S'i] Tr[5 4 527a] + Tr[S , 3 5i7 a ] Tr[5 , 47 A1 5 2 ] 
+ Tr[S 3 S'i] Tr[5 4 7^27a])GL(^i*) G^(fe3i.)* 

3 

+ 9 ~ £ f L {^(*3i *) ^2(fc3i *) (Tr[S 3 5i 7a ] Tr[5 4 5 2 ] 

+ Tr^] Tr[5 4 5 27a ]) G CT (fc 3 i*) G^(fc 31 *)* 
+^(A*l.) ^(fc 3 i*)(Tr[5 3 7a5i] Tr[5 4 5 2 ] 

+ Tr[S 3 S , i] Tv[S 4 S 2la })G a (k 31 ,)* G«Jk 31 *)} 

3 

+^ fahi *) J#(k3i *) (Tr[ 1 S 3 7 M 5 l7a ] TY[S 4 7^ 2 ] 

+ Tr[5 3 7 At 5' 1 ] Tr[547^27a])G^(fc 3 i*) G^(fc 31 *)* 

+ J>(A;31*) ^w( fe 31*) (Tr[<S'37Q!<S'l7nM] Tr^^T/d 

+ Tr[S 3 S l7u ]Tr[S 47a S 2llx }) G£"(*3i *)* G^(fc 31 *)} 

2 

{^(% *) ^,(*3i *) ^(hi *) (Tr[537 M S l7 5] Tr[5 4 5 2 7 5 ] 
+ Tr[S 3 5 l75 ] Tr[5 4 7 M S 2 75])G^(fc 3 i*) G n (k 3 i*)* 

+J r a(hi*) Fuifoi*) Flihi*) (Tr[S , 3 7 5 S'i7 At ] Ti[S^S 2 ] 
+ Tr[S 3 75^i] Tr[5 , 47 5 5 2 7 M ])G^(A ;31 *)* G^ 3 i*)} 
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In this formula, the J-(k) are the form factors given by Eq. (|T^) and k^i is the quadrivector for 
momentum transfer in the direct diagram k^ lJf = p^ — p^ 

The Mandelstam variables are related to the incident p\ * , P2 * and outcoming p% * , p^ * impul- 
sions by 



s 
u 



(Pi* +P2*) 2 
Ol* -P3*) 2 
G?l* -P4*) 2 



with s + t + it = 4 m 2 



m 



When the incident and outcoming particles are on their effective mass-shell pf 
take into account the energy momentum conservation p^ + = P3* +£>4*> we have 

s = 2 (ml +p!*.p 2 *) = 4 E 2 

t = 2 (m^-pi*.p 3 *) = 2 (ml -E 2 )(l- cos 6) 

u = 2 (ml-p^.p^) =2 (ml- E 2 )(l + cos 9) 



(63) 
(64) 
(65) 

2 and if 



(66) 
(67) 
(68) 



If the relative velocity between the background nuclear matter and the center of mass of the 
collision = (7, 71; sin a cos 92, 7U sin a cos 92,71; cos a) is finite, and when calculating terms 
containing the longitudinal part of the vector meson propagator, we will also have components 
along the vectors 



2 %i 
<?31 



and 



2 ?41 



We note the projections of the particle momenta on these directions pf = Pi-1]q^, etc. They 
are related to the Mandelstam variables defined above and to the components of the relative 
velocity through 



py 

pf 
pf 

P2 



P3 



Pf 



pX 



Ps 



u 



UJa 



+ "fv cos a 



jv COS OL\ 



u 



u — to, 



e 



s 


— m 2 






s 


— m 2 








— m 2 








— ml 







2 _ 
2 _ 
2 

Wffl 



The ojq and were defined in section (3.2). In terms of the Mandelstam variables, 

t 



LOq 



TO 4 



mi 



sin a sin (p\ / 1 



u 



cos a(l 



mi 



sin a sin <p\ 1 



t 



cos ail + 



-) 
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We split for convenience here the propagator of the to vector mesons in a transverse part G^t 
and a part along the vector 77: G U N(t) = G^iit) — G^rit). 

The spin-averaged scattering matrix for the Coulomb substracted elastic nucleon-nucleon scat- 
tering is given by: 



d(Tpp 

In 



dCnn 

dn 
1 

64ir 2 .t 



(69) 



with (e, d and i are the direct, exchange and interference terms) 



$4(7 = 

4(t — 

ci _ 
- 



gtJ*{t) \(4ml-t) 2 } G a (t)G* a (t) 



G a {u) G*(u) 



t 2 + u 2 - s 2 



+ 4m 2 (s- m 2 ) 



.{G a (t) G» + G*(t) G a (u)} 



4w 



>4w 



s d 



&4u) 



4 T~ 



(t) 4 (4 [(s - u) 2 + t(t + 8ml)] ■ {GMt)G*Mt)} 

+ 4 [At(p 2 e + pl e ) + 16p 2 e p 2 e + t 2 ] {G uN (t)G* U N(t)} 



+ 4 [4(s - u)piep 4 e + 4t(p 2 e + pl e ) + t 2 
x{G wJV (t)G^ r (t) + G: JV (t)G wr (t)}; 
<?^(u) 4 (4 [(a - t) 2 + «(« + 8m 2 )] . {GMv)Gl T {u)} 



+ 4 
+ 4 



MPi® +pI®) + 16 Pi®P4® + {G w jv(it)G£jv(«)} 



4(s - t)piep 4 ® + 4u(Pi® + #i®) 



x {G wiV (u)G* T (u) + Gl^G^u)} 



9i^l(u)Tl(t){- 8(s-2ml)(6ml-s) . {G^G^u) + G^G^u)} 



8(4m 2 - s)(p le p 4e +Pi (B P3( B ) - 4t(p 2 e +p 2 4e ) - 4u(p 2 e +p§ e ) 

+ 32piep 4 ®Pi®P3® - 2ut(2r? e .r? e - 1) x 
x {G^^G^H + G^(t)G^H} 

- [4s(pi® -P3®) 2 - 8u(p? e + £>!©) + 4t(pi e + f> 3 ®) 2 - 8um 2 

- 16m 2 (p 2 e +p 2 ffi - 4 Pl ®p 3 ®) - 2(a - 2m 2 ) (s - u) + 2(2m 2 - t)(4m 2 - t) 
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+ 2llt(l + -p 3 g 



{G w jv(«)G£t(*) + Gl N (u)GMt)} 



4s(pi e -P4e) 2 - 8t(Pi e + pIq) + M^ie + P^ 



8tmi 



- Wml(pj e +p 2 4e - 4piep 4e ) - 2(a - 2m2)(s - t) + 2(2m^ - u)(4ml - u) 



+ 2ut(l + -( Ple - P4e f) 



>2a2u 



2ct2w 



9ig^(t)^(t) {-16mi(s - u)]. {G a (t)Gl T (t) + G*(t)G uT (t)} 



+ 



+ 



-64m^ leP4e . {G a (t)G* uN (t) + G;(i)G wJV (i)} 



16(4m2 - t)(- +p 2 1Q +pl Q ) + 128m2p leP4e 



.{G x (t)G* x (t)} 



glglTl(u)Tl(u) ([-16mfr - t)]. { G CT (it) G^, ^ (it) + G»G wT (iz)} 
+ [-64m2pi e p3e]. {G^G*^) + G*(u)G uN (u)} 



+ 



16(4ml - u)(| +p? e + 128m^i©p 3 © 



.{G x (iz)G* x (iz)} 



glglTl(t)Tl(u) ([4(24m^ _ 4m 2 fl _ i 0m 2 t + {G(T ( t)G * T(n) + G ;(t) G ^ T ( u )} 
+ 2u(s + u) + 4n(pi e + pg e ) + 8spi©p 3ffi x 
x{G CT (t)G^(n) + G;(t)G^(n)}) 
+S?fl£^(«)^(t) ([4(24m2 - 4rr& - 10m^ + u 2 )] . {G ff (u)G* uT (t) + G^G^t)} 
+ [2t(a + t) + 4t(p? e + p2 ) + 8s Ple p ie x 
x {G ff (u)G^(t) + G£(u)G wJV (t)} 



9A^(^<7(«)^(tt)^(t) [32m,, (pi e + p 4e )(pi© + £>3©) + 8ut r/ e .r/ a 
x {G x (n)G* x (t) + G*(n)G x (t)} 



g ffS it^(t)^(t) ([32m, (( s - u ) Ple + t P4e )] . {G x (i)G* T (t) + G* x (t)G wT (t)} 

+ [16m* (pie +p 4e )(t + 4 me p 4e )] ■ {G x (t)G^(t) + G* x (t)G wJV (t)}) 

g ff g&F ff (u)J#(u) ([32m, ((a - t)p le + up 3 ©)] . |G x ( U )G: t ( U ) + G* x (u)G wT (u)} 

+ [16m* (pi© +p 3ffi )(n + 4pi e p 3ffi )] . {G x (n)G* ^(n) + G* (it)G w jv(u)}^ 
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-gaglT^F^Tliu) ([l6m*( Ple +pas)(u - t + 2ml)] ■ {G x (t)G* T H + G* x (t)G wT («)} 
+ [8m*(pi e +P4e)(4pi©p 3 e +«)] • {G x (i)G* ^(u) + G* x (t)G wN {u)} 



g a glF a {u)^(u)Fl(t) ([l6m„(pi© +Pse)(i - ^ + 2mji)] . {G x (n)G: T (t) + G*(u)G wT (t)} 



+ 



[8m* (pi® +p 3 e)(4pieP4e + *)] • {G x (n)G* Ar (t) + G* x (u)G uJV (t)}) 



'Scrlu "Stria; 



-&gjr*(t)F u (t) [l6m*(4m 2 - *)(Pie +P4e)] • {G x (t)G£(t) + G* x (t)G ff (t)} 
- 5 ^3 (n) ^ (n) [ 16m#(4m 2 _ n )( piffi +p3ffi )] . {G x (u)G» + G*»G CT (u)} 

(^(t)Jv(«)^ w («) [8m*(4m 2 - t)(p le +p 3e )] . {G CT (i)G* (u) + G*(t)G x (u)} 
+ J*{u)F a {t)F u {t) [8m*(4m 2 - u)( PlQ +p 4Q )] . {G a (u)G* x (t) + G* a (u)G x (t)}) 



gtT${t)[t 2 \Gn{t) G w (t)* 
^^ 4 (n) [u 2 ] G w («) G w (u)* 



4 ^(t) ^(u) [-].{G v (t) G v (ur + G v (ty G n (u)} 



$2n2a + S2-k2ct + ^27r2<T 






gl gl (j*(t) Fl(u) [ (4 ml ~ t)U ].{G a (t) G 7r (u)* + G a (i)* G»} 



+ ^(u)^(t) [ 



2 

(4 m 2 - u)t 



{G CT (u) G w (t)* + G ff (u)* G w (t)} 



^27r2w + ^27r2w + ^27r2w 
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5? 



2-k2w 



-glglFl{t)Fl{u)([Au{2ml - u)] . {G^G* (u) + Gl T {t)G^(u)} 
+ [8(4ml - s)p leP4e - 4t(p? e + pl e ) + 2ut x 
x {G wN (t)Gl(u) + G^G^u)}) 
-gl 9 lTl{u)Tl{t) {[4t(2ml -t)] . {G U T(u)G*(i) + G* T (u)G^(t)} 



+ 



8(4m^ - s)pi©p t © - 4n(p? e + p§ ffi ) + 2ut 
x {G wA r( U )G;(t) + G^^G^t)}) 



•Sfirlerlu; = 

(f a {t)^(t)^(u) [8m*u(p le +p 4e )] . {G x (t)G;H + G* (t)G w («)} 
Jv(«)^ w (u)^(t) [8m*t(pi e +Ps©)] • {G x (u)G;(i) + G* (u)G w (t)}) 
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Figure captions 

Fig. 1: N-N scattering through exchange of one dressed n, a, u meson 

Fig. 2: Fitting the experimental total elastic pp cross section: comparison of the performance 
of the various parameter sets. Experimental data points are indicated by the diamonds. The 
data set A, which is used for all calculations in this paper, corresponds to the full line. 

Fig. 3: Fitting the experimental differential elastic cross section: (a) comes from a smoothed 
version of the compilation of existing experimental data for five values of the beam energy 
T lab = 100, 300, 500, 700 and 900 MeV. The dashed part on the left hand side of the figure 
corresponds to the range where Coulomb scattering is dominant and which was not taken into 
account for fitting, (b) Result of the fit for parameter set A and five values of the beam energy. 
(c) Differential cross section at saturation density and vanishing temperature, for five values of 
the beam energy. 

Fig. 4: Total elastic p-p cross section as a function of the beam energy at T=0, several 
densities. From top to bottom: n/n sat = 0, 0.05, 0.1, 0.2, 0.5, 1., 2., 4. 
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Fig. 5: Ratio of the value of the in-medium cross section to its value in the vacuum as a 
function of density n/n sa t (log. scale) for seven values of the temperature T=0, 40, 80, 120, 
160, 200 and 240 MeV. Cmed/^free is shown on (5a) for a beam energy Ti ab = 100 MeV and on 
(5b) for T lab = 300 MeV. 

Fig. 6: (a) Ratio of the value of the in-medium cross section to its value in the vacuum as a 
function of temperature for five values of the density n/n sa t=0., 0.1, 1., 2. and 4. and for two 
values of the beam energy Ti ab =100 MeV (full line) and T iab =300 MeV (dashed line), (b), The 
effective mass of the nucleon as a function of temperature is shown for comparison. 

Fig. 7: Ratio of the value of the in-medium cross section to its free value: Contour plot in the 
density-temperature logio(n/n sa t) vs. T plane, for three values of the energy of the incident 
particles T lab =100, 300 and 500 MeV. 

Fig. 8: Differential elastic p-p cross section at T = as a function of the scattering angle, for 
a value of the beam energy TJ ab = 800 MeV, and several densities: n/n sat = 0, 0.05, 0.2, 0.5, 
1., 2., 4. 

Fig. 9: Ratio of the "effective masses" (as defined in Eq. 27) to their value in the vacuum, 
and mixing parameter (Eq. 29), as a function of momentum transfer at T=0, n/n sat = 0.1. 
There is a critical value of the momentum transfer q = 2pp. Also shown are the values of the 
dynamical masses obtained by solving the dispersion relation at zero 3- momentum k = (u>f l , 0): 
tf" = uf, D(k) = 0. 

Fig. 10: Comparison of experimental constraints on the reduction factor on the in-medium 
cross section from BUU and AMD, Brueckner calculations, and our model. The BUU calcu- 
lation is represented by the calculation of Klakow and is indicated by two arrows for densities 
n/n sa t=l. and 0.5. The AMD calculation is represented by the results of a simulation by 
Tanaka et al. of p+ 9 Be (diamonds) and p+ 40 Ca (triangles) collisions, thought to probe den- 
sities n/n sat ~ 0.1 and 1. respectively. Theoretical results calculating the medium effects in 
the Brueckner approach are taken from the paper of Li and Machleidt at densities n/n sat = 
0.5 (dotted line) and 1. (full line) and vanishing temperature. The results of this work are 
displayed for two values of the temperature T = (thin lines) and T = 50 MeV (thick lines) 
and two values of the density n/n sat =0.1 (dashed lines) and n/n sa t =1- (dash-dotted line). 

Fig. 11: The uj-q plane, particle/hole damping region and zero sound branches. The zero- 
sound branches are calculated for g a = 3.8, = 9.3 without form factors. For a fixed value of 
the beam energy, scattering angle and quadri-momentum transfer, and by varying the relative 
velocity existing between the rest system of the fluid and the center of mass of the collision, 
the point of coordinates (uj-q) describes a line which may twice cross the zero sound branches. 
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